Impurity induced density of states and residual transport in nonunitary 

superconductors 



O 

o 

(N 

> 
O 

:z; 

(N 



O 

o 
I 

^ : 



O 

o 



> 

p 

od 
O 

o 



X 



T. R. Abu Alrub and S. H. Curnoe 

Department of Physics and Physical Oceanography, 
Memorial University of Newfoundland, St. John's, NL, A IB 3X7, Canada 

We obtain general expressions for the residual density of states, electrical conductivity and thermal 
conductivity for non-unitary superconductors due to impurity scattering. We apply the results to 
the so-called 'B phase' of PrOs4Sbi2, which we describe using a non-unitary gap function derived 
from symmetry considerations. The conductivity tensor has inequivalent diagonal components due 
to off-axis nodal positions which may be detectable in experiments. 
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I. INTRODUCTION 

Non-unitary pairing in superfluids was first described 
by Leggettfi but the Ai phase of '^He is the only well- 
established example of this, so far. However, recently 
non-unitary pairing was observed in the heavy fermion 
superconductor PrOs4Sbi2 by Aoki et al^ A physically 
significant consequence of non-unitary pairing is a lifting 
of the degeneracy of the superconducting energy gap, so 
that two different energy gap branches, both of which 
are anisotropic, are observable. Multi-gap behaviour has 
been observed in PrOs4Sbi2^i^i^i^ii2ii^ but so far this has 
mainly been attributed to multi-band superconductivity, 
and gap splitting due to non-unitary pairing has received 
little consideration, in spite of numerous citations of Aoki 
et aZ.'s results. 

Superconductivity in PrOs4Sbi2 is believed to be 
unconventional,^iHiiia5ii6a2ii8ii9i20 j,^^ ^^^^.^^ 

trons are in a spin triplet configuration^i^ and the super- 
conducting state has broken time reversal symmetry and 
is non-unitary.^ Low temperature power law behaviour, 
indicative of the presence of nodes in the gap func- 
tion, has been observed in thermodynamic and transport 
measurementS j^'^d^d^d^i^o but some experiments have 
found the gap function to be nodelessi^i^iiS^iS^ Other ex- 
periments observed two superconducting phases, possibly 
with different symmetries , ^d^^d^, 24^25, 26^27, 28, 29 g^ggegi;. 

ing a multi-component superconducting order parame- 
ter. These two phases are known as the "A phase" and 
the "B phase". If it exists, the A phase occupies only 
a small region of the phase diagram just below Hc2{T). 
Thus, most measurements, including those cited above, 
have probed the B phase. 

The three dimensional representation of the 
point group Th best describes superconductivity in 
PrOs4Sbi2i^i^ This representation yields several su- 
perconducting phases, of which four are accessible 
from the normal state by a second order phase tran- 
sition. We have previously identified the states 
02(02) X /C and D2{E), with order parameter compo- 
nents (0,0, Iryil) and (0,«j772|, I771I), as the A phase and 
B phase, respectively 1^1^ Here, 02(02) is the symme- 
try group with elements {E, Cf , C/(7r)C|, C/(7r)C|} while 



D2{E) = {£;,Cf/C,C/i(7r)C|/C,i7i(7r)C|}.3-a The corre- 
sponding gap functions are unitary for the A phase, with 
two point nodes in the [00±1] directions, and nonunitary 
for the B phase, with four nodes on unusual points on 
the Fermi surface, [0,±q;,±/3]. 



Low temperature transport is an effective probe for 
the symmetry of the gap functionJ^iiii^^i^i^i Impu- 
rities induce and scatter quasiparticles at the nodes 
and the conductance remains finite even in the limit 
of zero frequency and temperature. Usually, two lim- 
iting cases of impurity scattering are considered, the 
"Born limit" (weak scattering) and the "unitary limit" 
(strong scattering). The unitary limit is associated with 
non-magnetic substitutions of magnetic ions in heavy 
fermion superconductorsi^i^^i^i^ The self-energy due 
to isotropic impurity scattering is obtained from the T- 
matrix,'^^''^^ I](fc,cj) = {ni/TrNn)T{k, k,uj), where Ui is 
the impurity concentration, iV„ is the density of states in 
the normal state, and the T-matrix is the self-consistent 
solution to T{lu) ^ V + VOo{lu)T{ll;), where V is the 
impurity potential, Oq{uj) = {l/iTNn)^]^0{k,tjj) and 
0{k,ijj) is the electronic Green's function in the su- 
perconducting state. The self-energy is then S(w) — 
(ni/7riV„)G'o(w)/[c^ — G'g(w)], where c is related to the 
phase shift, c = cot 5q. In unitary limit c — > O, while c — > 
00 in the Born limit. The main result of this approach is 
a renormalisation of the frequency w — > w — iT[ijj) due to 
impurity scattering. We will use this result to find im- 
purity induced residual density of states and transport 
coefficients. 



The outline of this paper is as follows: in Sec. II we 
define the gap function, the mean field Green's functions 
and spectral functions. In Sees. HI, IV and V we de- 
rive general expressions for the impurity induced quasi- 
particle density of states, the electrical conductivity and 
the thermal conductivity in a nonunitary superconduct- 
ing state. In Sec. VI we apply our results to the nonuni- 
tary B phase in PrOs4Sbi2, and we summarise our results 
in Sec. VII. 
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II. MEAN FIELD RESULTS 



where 



In the following we state the main results of the mean 
field treatment of an effective pairing Hamiltonian (see 
Ref.lUfor details). 

The gap function is a 2 x 2 matrix in pseudospin space. 
For triplet pairing it can be parametrised in terms of an 
odd pseudovectorial function d{k) as 



Afc = i[cr-dk](7y = ( 



+ idy{k) dz{k) 
(fc) dx{k) + idy{k) 



(1) 

When AfcA^ is proportional to the unit matrix, the pair- 
ing is said to be "unitary" . Non-unitary pairing occurs 
only in the triplet channel and only when = idk xd^ ^ 
0. Non-unitary states necessarily have broken time rever- 
sal symmetry. However, note that, for example, pairing 
of the form d^ = (fcx + iky)z (proposed for Sr2Ru04) 
breaks time reversal symmetry but is unitary. The quasi- 
particle energies are 



Afc± = [\dk\^ ± Igfel 



,1/2 



(3) 



Thus, non-unitary pairing lifts the gap degeneracy. 



For triplet pairing, the normal and anomalous quasi- 
particle Green's functions are^MS 



G{k,iujn) 



iujn + £k\ 



(4) 



■ ^ [ujI + el + \dk\'^]dk - iqk X dk — -, 



k± 



[4 



1/2 



(2) It is useful to expand these expressions as 
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where 



Uk±Vk± 



£fc 1 2 1 ri ^'^ 1 

^^^ = 2[1-^] 



2-Bfc±^ 



(8) 



are the extended coherence factors for this particular 
state. Note that the following identity has been used 
in deriving the above expressions 



be included by replacing iw„ with iw„ — E(zw„). The 
retarded self-energy is Sret(w) = S(iw„ ^ cu + i6) = 
—iT{uj) where the real part is assumed to be frequency 
independent and absorbed in the chemical potential. 

The spectral function A'^{k,uj) (and similarly 
A^{k,uj)) is defined by 



G[k, luon) = I duj ■ 



10J„ 



(10) 



i{qk xdk)-(T = [qk ■ (T)(dk cr) -qk -dk 

= {Qk ■ o-)(dk ■ a-) (9) 

where qk ■ dk ~ because qk -L dfc. The self-energy can 



Usually, the spectral function is just — i3G"'^*(A;, cu), but 
in this case, because the Green's function has a complex 
numerator, the spectral function must be extracted more 
carefully. Using ^ and ([7]), one finds 
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27r 
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Ak+ [(Lu-Ek+r+T^iLo) iL^ + Ek+r+T^Lo)_ 



(12) 



with the spectral functions in hand, we can proceed to 
calculate the density of states and the transport coeffi- 
cients. 



III. DENSITY OF STATES 

The quasiparticles density of states can be defined in 
terms of the spectral function as 



iVM = ^Tr[A«(fc,a;)] 



(13) 



using pT|) we find the general expression for the density 
of states in a nonunitary superconductor, 

N{lo) = H±^(^ -Ek±) + vI±6(lo + Ek±)] (14) 

k,± 

in the absence of impurities. For small lu, in the vicinity 
of the gap node, we have Vk± ~ 0, Uk± ~ 1, and (|14p is 
reduced to^ 



7V(l^) «^5(c^-£;fe±). 



(15) 



fc,± 



When the impurities arc included the density of states 
becomes 



TT ^ ^ 



k,± 



1 



{Lu-Ek±r + r^Lu) 



It is clear from (|T6l) that the residual density of states 
depends on the impurity concentration through the self- 
energy r(o). 



IV. ELECTRICAL CONDUCTIVITY 



The DC electrical conductivity is defined by the Kubo 
formulai^ 



lim 



(17) 



where 

n(q, iil,, 



dre^'^-^ (r.jg(r)j-„(0)) (18) 



is the current-current correlation function. The electrical 
current is defined by 



h.s 



4+,,.(^)cfe,.(r). (19) 



The current-current correlation function is therefore 

k 

+F{k, iuJn)F\k + q, ztj„ -f ^^]„)]• (20) 

The conductivity vanishes when the self-energy is absent, 
and the contribution from the anomalous part vanishes 
even when the sclf-encrgy is included. In the limit q — > 
the correlation function is 



(16) n(il7„) = e2 E VfVf), E Tr[G(fe, ltOn)G{k, ItUn + l^n)] 



_ ^ (21) 

To evaluate this correlation function we follow the ap- 
proach of Refs. [U and fill and rewrite the Green's func- 
tion in terms of the spectral function (jlip and sum over 
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Matsubara frequencies. This eventually leads to 

/oo 
du)' 

[npiLo') -npioj' + ^)]. (22) 
Then the DC electrical conductivity (fTTl) is 

/oo 
rfcj'Tr[A^(c^')^fe(^')] 
-oo 

(23) 

= 5{uj'), and the 



k 



duo' 



In the limit T ^ we have - ^"J^^ = 
conductivity is 

5 = 7re2^t;;^t;;^Tr[I^(o)IG(o)]. 



(24) 



Using (fTTj) we finally obtain the conductivity for a non- 
unitary superconductor, 



~ ^ ^ \ 

(7 = > VpVp 

n ^-^ 

k 

where To r(tj 0). 



(25) 



As in the electrical conductivity, the anomalous part does 
not contribute to the thermal conductivity. Finally, the 
correlation function is expressed in terms of the spectral 
function as 



3n,et(i^) 



k 

Tr[I«(c.' + f^)I«(-c^')] 
+ 2 



(29) 



Substituting this into the Kubo formula ((26l) and evalu- 
ating in the limit f2 — + and T — > 0, wc find 



T 



(30) 



Comparing ([^1]) and ([50]) we can see that the 

. 9 

is satisfied. Explicitly, 



Wiedemann-Franz law -% - 
the thermal conductivity is 



3e2 



T 3 



k 



(31) 



VI. APPLICATION TO PrOs4Sbi2 



V. THERMAL CONDUCTIVITY 

The DC thermal conductivity is defined by the Kubo 
formulaiii 



K 1 ^Ilreti^) 

— = ^ lim . 

T T2 n^o Q 



(26) 



The heat current can be written in second quantization 
form as 



Mr) 



2m* ^ 

k,s 



t[k + q] — Ck+q,s{T) 



(27) 



This form is similar to (4.17) in Rcf. [Tl| except that we 
have neglected the term proportional to the gap veloc- 
ity, which we assume to be much smaller than the Fermi 
velocity. The current-current correlation function is then 



man) = ^ E 



VpVp 



Tr[G'fc(iw„ -I- iQn)Gk{-iuJn) 
~Fk{-iuJu)Fl{ioJn + ^^^„)] (28) 



As discussed in the Introduction, we assume that the 
gap function for the A phase is 



^k^\m\[aX + ^"kf^ 



(32) 



where a and h are undetermined constants, while for the 
B phase it has the form 



A 



k± 



±2 



I V'^ 

\Vi\\m\\kx\^a-^b-^kl+a'^kl + h^kl . (33) 



which is non-degenerate The gap function in the 
A phase is unitary and has two cusp point nodes in the 
[00±1] directions. The lower branch of the B phase gap 
function has four point nodes which are in the fcj, = 

plane at the positions \/\rii\'^b'^ — \il2\'^ci'^kx — ±\r]2\bkz if 
?7ip62 |j^^|2q2. gjgg ^j^gy g^j.g jj-^ kz = plane. We 

will assume the former in our calculations. Since we are 
interested in the very low temperature regime, we will 
consider only the B phase. 

The gap function of the B phase in the vicinity of nodes 
can be linearised as 



1.12 



(34) 
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where v = v^|?7ip6^ — |?72pa^, K = ^ ky and 



\Vi\b 



\m\a 
\m\h 



(35) 



/c|| and k±_ (used below) are momenta parallel and per- 
pendicular to the Fermi surface at the node. The upper 
branch, which is degenerate with the lower branch on the 
line kx ^ Q between each pair of nodes, is properly in- 
cluded with this linearisation of the gap function. There- 
fore, we relabel the two branches of the gap function as 
shown in Fig. 1. Thus for any function we have 



f{E+) + f{E^)^f{E^) + f{E2). 



(36) 



Each branch 1 and 2 has two cusp point nodes and the 
contribution to the excitation spectrum from each branch 
is equal. With this picture in mind, we now calculate the 
density of states and the transport coefRcients. 





FIG. 1: (Color online) Gap function for the B phase of 
PrOs4Sbi2 drawn in the k^-k^ plane over a spherical Fermi 
surface (bold black). Left: the branch is shown in blue 
(dashed) and the '-' branch in red (solid). Right: the '1' 
branch is shown in blue (dashed) and the '2' branch in red 
(solid). 



A. Density of states 

The density of states was calculated previously in Ref. 
I31I in the absence of impurities; here wc will include the 
effect of impurities starting from (|16p . Linearising the 
gap function as described above, we find 



iV(.) = ^M2^ 

TT ^ ^ 



(Pk 



1 



(27r)3 {oj-Ekf+T^{oj) 



(37) 



where there is a factor of 2 because there are two branches 
of the gap function and the sum is over the two nodes 
in each branch. To perform the integration we change 
variables to = w^(fcj^ -I- ky) -f v']^k\ 

2F(w) h 1 

(38) 



N{uj) 



El 



Pa 



dp 



and introduce a cutoff po- Performing the integration we 
arrive at the result 



1 



1 ( Po - ^ 



tan 

-\-LoT{uj) In 
+Por(w) 



tan 



[po-iuf+T^{uj) 
u;2+r2(a;) 



Setting r(a;) = we obtain our previous resuHjai 



_ b 2w^ 



(39) 



(40) 



which has a quadratic dependence on frequency as ex- 
pected for point nodes. In the limit w ^ (|39p becomes 



N{0) = 



2 b 



^ 



tan 



To 



To 

... 

This is the zero energy density of states induced by impu- 
rities. The cut-off is normally taken to be the size of the 
Brillouin zone^^ but it may be more physical to use the 
reciprocal of the range of the single impurity potential, — 
Po oc A^^. In terms of the ratio (po/Tq) the two limits 
are 



— 3> 1 (unitary) 
To 
Po_ 

To 



< 1 (Born) 



In the unitary limit the density of states is 
2 b poT^ 



iV(0) 



[\m\^b^ - \m\^a^]vF 



(42) 
(43) 

(44) 



where u refers to unitary scattering. If Fc is the criti- 
cal scattering rate at which the superconductor becomes 
normal at the node, then we can write (|44p as 



^(0) 



imp 



(45) 



In the Born limit, the density of states vanishes as Fq. 

The presence of residual density of states, in general, 
gives a contribution linear in temperature to the specific 
heat and the nuclear spin relaxation rate at low tem- 
perature. The prefactor dependence on impurity doping 
may be helpful in identifying the symmetry of the order 
parameter. The specific heat is^ 



C(T) = 



T 



duJUJ^N{uj) 



duj 



At low temperature this yields 

{C{T)/T) f; 



{C{T)/T\ 







(46) 



(47) 
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and the nuclear spin relaxation rate if 



,39 



IT 



{1/TTi)t 

(l/TTi)„ 



dujN{uj)N{uj -ujo) 



doj 



(48) 



Electrical and thermal conductivities 



Beginning with and making use of ([M]) , we divide 
the integration into four parts, each centred about one 
node in the gap function. The factor vpVp is evaluated 
at each node; the sum over nodes yields 



a^x and azz due to the off-axis nodal positions and the 
choice of a particular domain of superconducting phase. 
This domain is represented by order parameter compo- 
nents (0, z|772|, |?7i|)- If all six domains are present then 
all diagonal components of the conductivity tensor will 
be equal. The a^x component is proportional to the 
parameter |?72| which is absent in the unitary A phase. 
Therefore, measurement of residual conductivities in a 
domain-pinned set-up, such as the one used in direc- 
tional dependent thermal conductivity measurementai^ 
could determine the direction of nodes. Of all the pos- 
sible SC states in tetrahedral systems, D2{E), with OP 
components (0, z|772|, |f7i|), is the only one with off-axis 
nodes.— 1^ 



4 

E 



vpVp — Avp 



WW " 






(49) 



The remaining integration is the same for each part. Per- 
forming the same change of variables as in the density of 
states calculation, we find 



Y.. = l '"fVf 



dp 



27r3 a[\T]i\%^ -\r,2\'^a^]vF Jo [p^ + Tlf 
and completing the integration we get 



(50) 



cr = — VfTq 



-1 



tan 

a|r)2 



(Po/ro) 



1 + (Po/Toy 



6|r,i|2[62|^^|2 






-a^ I j;2 P] 







(51) 



This is the impurity induced DC electrical conductivity 
for the B phase of PrOs4Sbi2. The thermal conductivity 
can be easily obtained by using the Wiedemann-Franz 
law. In the unitary limit (^ 3> 1), the term which in- 
cludes tan^^ (1%) ~ ^ ^^^^ dominate, the conductivities 
become 



and 



_ pit 








/ /'l'?2|'' 

b\ni\'''[b'-'\rn\''-o.-'M'^] 





(52) 



V 








1 



(53) 

Thus the conductivities in the B phase of PrOs4Sbi2 
are non-universal (dependent on impurity concentration) 
for unitary scattering but vanish in the Born limit. The 
conductivity tensor has two inequivalent components. 



C. Discussion 

There have been several studies on Ru and La doped 
samples (Sii^ii^'^i^ with the surprising result that Ru 
substitution leads to a doping-dependent residual den- 
sity of states and resistivity , -'^^''^^ while La substitution 
does noti^ In PrOs4Sbi2, it is speculated that quadrupo- 
lar fluctuations of the Pr ions play a role similar to the 
magnetic fluctuations of Ce and U ions in other heavy 
fermion superconductors, thus substitution of the Pr ions 
by La would be expected to produce unitary scatterers. 
However, in contrast to Eq. [48l there is no dependence 
on doping on NQR relaxation rate beyond the La con- 
centration X — 0.05. 

Both Pri„j;Laa;Os4Sbi2 and Pr(Osi_a;Ru2;)4Sbi2 are 
superconducting for the entire range of a;, and both 
become s-wave superconductors at some intermediate 
value of X. In the Ru doped series, Tc has a mini- 
mum at a; = 0.6, with a leveling off of the specific heat 
at the same value. This suggests that a phase transi- 
tion between triplet and singlet superconductivity oc- 
curs at X ~ 0.6, with possibly a region of co-existence 
of these two phases.'^ A 0.4% change in lattice constant 
occurs between PrOs4Sbi2 and PrRu4Sbi2,''^ and effects 
due to quadrupolar fluctuations appear to be absent in 
PrRu4Sbi2. In the La doped series, Tc decreases linearly 
along the entire range of x, while the specific heat levels 
off at x « 0.3. 

According to and ((5^ . the dependence of the 

residual density of states and resistivity on Ru doping 
suggests that the scattering from Ru ions is unitary. Uni- 
tary scattering due to the substitution of Os by Ru may 
be explained by noting that quadrupolar fluctuations of 
the Pr ions are charge density fluctuations and will cou- 
ple to, and possibly be enhanced by, quadrupolar lattice 
vibration modes. The change in lattice constant that ac- 
companies Ru doping will alter the quadrupole moment 
of those modes. In addition, Ru substitution has a strong 
effect on the low-lying crystal electric field (CEF) lev- 
els of the Pr ions which eventually removes quadrupole 
fluctuations!^ La substitution produces a much smaller 
change in lattice constant and has a much weaker effect 
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on the Pr CEF levels. Nevertheless, it is still difficult 
to explain why there is no dependence at all on the La 
concentration. 

VII. SUMMARY AND CONCLUSIONS 

It is evident from (fTHl) . ([25|) and ([3T|) that the main ef- 
fect of a non-unitary superconducting state is a lifting of 
the gap degeneracy, and that this would be observed as 
multi-gap behaviour similar to what could be expected 
for multi-band superconductivity. There are, however, 
some differences which we outline here. We base the fol- 
lowing discussion on the unitary state 02(02) x IC and 
the non- unitary state D2{E), with order parameter com- 
ponents (0,0, |?7i|) and (0, z|772|, respectively. There 
are many other states, but all the rest are either nodeless, 
or else they have a C3 symmetry element which has been 
positively ruled out by experiment .-iS. 

In a multi-band superconductor with a single Tc the 
symmetry of the superconducting order parameter should 
either be the same on both bands, or possibly, super- 
conductivity on one band is a secondary order parame- 
ter to superconductivity on the other. The alternative, 
which is the simultaneous appearance of two different 
order parameters, would be unprecedented. This means 
that the symmetries of superconducting states on the dif- 
ferent bands should either be the same, or have a group- 



subgroup relation. For example, in MgB2, the archetypal 
multi-band superconductor, s-wave superconductivity is 
observed as a full gap for both bands. The best can- 
didates for nodal superconductivity in the triplet chan- 
nel in PrOs4Sbi2 are the unitary state 02(02) x K, and 
the non- unitary state 02(E), and neither of these has 
secondary order parameters.— Therefore multi-band su- 
perconductivity entails nodes at the same places for both 
gaps, unless that part of the Fermi surface is missing. On 
the other hand, the non-unitary superconducting state 
has nodes in the lower branch and a fully gapped upper 
branch. This difference may help to distinguish these two 
possibilities. 

To summarise, we have found general expressions for 
the residual density of states and electrical and ther- 
mal conductivities due to impurity scattering, and we 
have applied the results to the non-unitary B phase of 
PrOs4Sbi2. The nodal positions of the non- unitary state 
02(E) are unique among all the superconducting states 
for crystals with tetrahedral symmetryj^SiSi in that they 
are not found on a symmetry axis. Inequivalent diago- 
nal components of the conductivity tensor would be an 
unmistakable signature of such a state. 

Acknowledgments 

This work was supported by NSERC of Canada. 



1 A. J. Leggett, Rev. Mod. Phys. 47, 331 (1975). 

^ Y. Aoki, A. Tsuchiya, T. Kanayama, S. R. Saha, H. Sug- 

awara, H. Sato, W. Higemoto, A. Koda, K. Ohishi, K. 

Nishiyama, and R. Kadono, Phys. Rev. Lett. 91, 067003 

(2003). 

^ D. E. MacLaughlin, L. Shu, R. H. Heffner, J. E. Sonier, F. 
D. Callaghan, G. D. Morris, O. O. Bernal, W. M. Yuhasz, 
N. A. Frederick and M. B. Maple, Proceedings of the Inter- 
national Conference on Strongly Correlated Electron Sys- 
tems, Houston, May 2007 (unpublished). 

* C. S. Turel, J. Y. T. Wei, W. M. Yuhasz and M. B. Maple, 
Physica C, 463-465, 32 (2007). 

^ G. Seyfarth, J. P. Brison, M.-A. Measson, D. Braithewaite, 
G. Lapertot and J. Flouquet, Phys. Rev. Lett. 97, 236403 
(2006). 

® M. Yogi, T. Nagai, Y. Imamura, H. Mukuda, Y. Kitaoka, 
D. Kikuchi, H. Sugawara, Y. Aoki, H. Sato and H. Harima, 
J. Phys. Soc. Jpn. 75, 124702 (2006). 

T. Sakakibara, A. Yamada, J. Custers, K. Yano, T. 
Tayama, H. Aoki, and K. Machida, J. Phys. Soc. Jpn. 76, 
051004 (2007). 

* M.-A. Measson, D. Braithewaite, J. Flouquet, G. Seyfarth, 
J. P. Brison, E. Lhotel, C. Paulsen, H. Sugawara and H. 
Sato, Phys. Rev. B 70, 064516 (2004). 

® G. Seyfarth, J. P. Brison, M.-A. Measson, J. Flouquet, K. 

Izawa, Y. Matsuda, H. Sugawara, and H. Sato, Phys. Rev. 

Lett. 95, 107004 (2005). 
° P. A. Lee, Phys. Rev. Lett. 71, 1887 (1993). 
^ A. C. Durst and P. A. Lee, Phys. Rev. B 62, 1270 (2000). 



M. B. Maple, E. D. Bauer, V. S. Zapf, E. J. Freeman, N. 

A. Frederick and R. P. Dickey, Acta. Phys. Pol. 32, 3291 
(2001). 

E. D. Bauer, N. A. Frederick, P.-C. Ho, V. S. Zapf, and M. 

B. Maple, Phys. Rev.B 65, 100506(R) (2002). 

K. Izawa, Y. Nakajima, J. Goryo, Y. Matsuda, S. Osaki, 
H. Sugawara, H. Sato, P. Thalmeier, and K. Maki, Phys. 
Rev. Lett. 90, 117001 (2003). 

E. E. M. Chia, M. B. Salamon, H. Sugawara, and H. Sato, 
Phys. Rev. Lett. 91, 247003 (2003). 

A. D. Huxley, M.-A. Measson, K. Izawa, C. D. Dewhurst, 
R. Cubitt, B. Grenier, H. Sugawara, J. Flouquet, Y. Mat- 
suda and H. Sato, Phys. Rev. Lett. 93, 187005 (2004). 
M. Nishiyama, T. Kato, H. Sugawara, D. Kikuchi, H. Sato, 
H. Harima, and G.-q. Zheng, J. Phys. Soc. Jpn. 74, 1938 
(2005). 

N. A. Frederick, T. A. Sayles and M. B. Maple, Phys. Rev. 
B 71, 064508 (2005). 

W. Higemoto, S. R. Saha, A. Koda, K. Ohishi, R. Kadono, 
Y. Aoki, H. Sugawara and H. Sato, Phys. Rev.B 75, 
020510(R) (2007). 

K. Katayama, S. Kawasaki, M. Nishiyama, H. Sugawara, 
D. Kikuchi, H. Sato, and G.-q. Zheng, J. Phys. Soc. Jpn. 
76, 023701 (2007). 

D. E. MacLaughlin, J. E. Sonier, R. H. Heffner, O. O. 
Bernal, B.-L. Young, M. S. Rose, G. D. Morris, E. D. 
Bauer, T. D. Do and M. B. Maple, Phys. Rev. Lett. 89, 
157001 (2002). 

H. Suderow, S. Viera, J. D. Strand, S. Bud'ko and P. C. 



8 



Canfield, Phys. Rev. B 69, 060504 (2004). 

H. Kotegawa, M. Yogi, Y. Imamura, Y. Kawasaki, G.-q. 

Zheng, Y. Kitaoka, S. Ohsaki, H. Sugawara, Y. Aoki, and 

H. Sato, Phys. Rev. Lett. 90, 027001 (2003). 

Y. Aoki, T. Namiki, S. Ohsaki, S. R. Saha, H. Sugawara 

and H. Sato, J. Phys. Soc. Jpn. 71, 2098 (2002). 

R. VoUmer, A. Fai/3t, C. Pfleiderer, H. v. Lohneysen, E. 

D. Bauer, P.-C. Ho, V. Zapf, and M. B. Maple, Phys. Rev. 

Lett. 90, 057001 (2003). 

T. Tayama, T. Sakakibara, H. Sugawara, Y. Aoki and H. 
Sato, J. Phys. Soc. Jpn. 72, 1516 (2003). 
P.-C. Ho, N. A. Frederick, V. S. Zapf, E. D. Bauer, T. D. 
Do, M. B. Maple, A. D. Christianson, and A. H. Lacerda, 
Phys. Rev. B 67, 180508(R) (2003). 

N. Oeschler, P. Gegenwart, F. Weickert, I. Zerec, P. 

Thalmeier, F. Steglich, E. D. Bauer, N. A. Frederick and 

M. B. Maple, Phys. Rev. B 69, 235108 (2004). 
=^ K. Grube, S. Drobnik, C. Pfleiderer, H. v. Lohneysen, E. D. 

Bauer and M. B. Maple, Phys. Rev. B 73, 104503 (2006). 
^° I. A. Sergienko and S. H. Curnoe, Phys. Rev. B, 70, 144522 

(2004). 

T. R. Abu Ahub and S. H. Curnoe, Phys. Rev. B 76, 
054514 (2007) 

P. J. Hirschfeld, W. O. Putikka, and D. J. Scalapino, Phys. 
Rev. Lett. 71, 3705 (1993); P. J. Hirschfeld, W. O. Putikka, 
and D. J. Scalapino, Phys. Rev. B 50, 10 250 (1994); P. J. 
Hirschfeld and W. O. Putikka, Phys. Rev. Lett. 77, 3909 
(1996). 



M. J. Graf, S-K. Yip, J. A. Sauls, and D. Rainer, Phys. 
Rev. B 53, 15 147 (1996). 

A. V. Balatsky, A. Rosengren, and B. L. Altshuler, Phys. 
Rev. Lett. 73, 720 (1994). 

C. J. Pethick and D. Pines, Phys. Rev. Lett. 57, 118 
(1986). 

B. Arfi and C. J. Pethick, Phys. Rev. B 38, 2312 (1987). 
S. Schmitt-Rink, K. Miyake, and C. M. Varma, Phys. Rev. 
Lett. 57, 2575 (1986). 

P. J. Hirschfeld, P. Wolfle and D. Einzel, Phys. Rev. B 37, 
83 (1988). 

M. Sigrist and K. Ueda, Rev. Mod. Phys. 63, 239 (1991). 
*° There is a minus sign in front of the flrst term in Q which 

is missing in Ref. [391 . 
'^^ G. D. Mahan, Many-Particle Physics (Plenum Press, 3rd 

edition, New York) (2000). 
■'^ N. A. Frederick, T. D. Do, P.-C. Ho, N. P. Butch, V. S. 

Zapf and M. B. Maple, Phys. Rev. B 69, 024523 (2004); 

N. A. Frederick, T. A. Sayles and M. B. Maple, Phys. Rev. 

B 71, 064508 (2005). 

C. R. Rotundu, P. Kumar and B. Andraka, Phys. Rev. B 
73, 014515 (2006). 

M. B. Maple, N. A. Frederick, P.-C. Ho, W. M. Yuhasz and 
T. Yanagisawa, J. Superconductivity and Novel Magnetism 
19, 299 (2006). 

L A. Sergienko, Phys. Rev. B 69, 174502 (2004). 



